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ABSTRACT 

A large number of cosmological studies now suggest that roughly two-thirds of the critical energy 
density of the Universe exists in a component with negative pressure. If the equation of state of such an 
energy component varies with time, it should in principle be possible to identify such a variation using 
cosmological probes over a wide range in redshift. Proper detection of any time variation, however, 
requires cosmological probes beyond the currently studied range in redshift of ~ 0.1 to 1. We extend our 
analysis to gravitational lensing statistics at high redshift and suggest that a reliable sample of lensed 
sources, out to a redshift of ~ 5, can be used to constrain the variation of the equation of state, provided 
that both the redshift distribution of lensed sources and the selection function involved with the lensed 
source discovery process are known. An exciting opportunity to catalog an adequate sample of lensed 
sources (quasars) to probe quintessence is now available with the ongoing Sloan Digital Sky Survey. 
Writing w(z) rs w + z(dw / dz) , we study the expected accuracy to which the equation of state today 
w Q and its rate of change (dw/dz) can simultaneously be constrained. Such a determination can rule 
out some missing-energy candidates, such as classes of quintessence models or a cosmological constant. 

Subject headings: cosmology: theory — cosmology: observations — gravitational lensing — quasars: 
general 



1. INTRODUCTION 

There is now strong evidence that the matter density of 
the Universe, both baryonic and dark, is smaller than the 
critical value predicted by inflation. Some of the same cos- 
mological probes that suggest the ratio of matter density 
to critical density is less than one prefer the presence of 
an additional component with energy density such that the 
total energy density of the Universe is in fact the critical 
value (e.g., Lineweaver 1998; Perlmutter et al. 1998; Riess 
et al. 1998). Type la supernovae and other observations 
have recently provided evidence that the equation of state, 
w, of this component is — 1 <; w = P/p <^ — 1/3 (Garnavich 
et al. 1998; Waga & Miceli 1998; Cooray 1999). 

The oldest known candidate for this additional energy 
is the cosmological constant, characterized by w = —1. 
However, other possibilities have now been considered, in- 
cluding a slowly-varying, spatially inhomogencous scalar 
field component better known as quintessence (e.g., Huey 
et al. 1998; Zlatev et al. 1998; Wang & Steinhardt 1998). 
Certain quintessence models, called tracker models (Fer- 
reira & Joyce 1997; Zlatev et al. 1998; see Steinhardt et 
al. 1998 for a comprehensive review), have the feature 
that the energy density of the scalar field, Q,q, is a fixed 
fraction of the energy density of the dominant component. 
Therefore such models may explain the coincidence prob- 
lem — why Q,q is of the same order of magnitude as £Im 
today. In these tracker models w is a function of redshift, 
and typically varies from w ~ 1/3 during the radiation 
dominated era to w ~ —0.2 during the matter-dominated 
era and finally to a value w ~ —0.8 during late epochs 
(today). Time variation of the equation of state is even 



more prominent in scalar field models involving pseudo- 
Nambu-Goldstone bosons (PNGB models) (Frieman et al. 
1995). In such models the field is frozen to its initial value 
due to the large expansion rate, but becomes dynamical 
at some later stage at redshift z x . Likely values for z x 
are roughly between 3 and (Coble et al. 1997), which 
means that interesting dynamics — and hence the vari- 
ation in the equation of state — happen at redshift of 
a few. Huterer & Turner (1999) point out that distance 
measurements to Type la supernovae offer a possibility to 
reconstruct the quintessence potential, while Starobinsky 
(1998) and Wang & Steinhardt (1998) suggested the possi- 
bility of using cluster abundances as a function of redshift 
as well. 

In § 2, we study constraints on w{z) based on current 
Type la supernovae distances, gravitational lensing statis- 
tics, and globular cluster ages. In § 3, we consider the 
possibility of imposing reliable constraints on w(z) based 
on cosmological probes at high redshift, in particular grav- 
itational lensing statistics. We follow the conventions that 
the Hubble constant, Ho, is 100 h km s _1 Mpc" 1 , the 
present matter energy density in units of the closure den- 
sity is fijif) and the normalized present day energy density 
in the unknown component is Qq. 



2. CURRENT CONSTRAINTS ON W(Z) 

Since a generic quintessence model has a time-varying 
equation of state, one should be able to distinguish it from 
models where the equation of state is time independent 
(e.g, Turner & White 1997). In this Letter, we write the 
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equation of state of the unknown component as 

w(z) w w + z(dw/dz) . (1) 

This relation should be a good approximation for most 
quintessence models out to redshift of a few and, of course, 
exact for models where w(z) is a constant or changing 
slowly. Note that negative (dw/dz) corresponds to an 
equation of state which is larger today compared to early 
epochs. Models where the scalar field is initially frozen 
typically exhibit such behavior, while for tracker field mod- 
els (dw/dz) > 0. 

In order to constrain w(z), we extend current published 
analyses which have so far considered the existence of a 
redshift-indepcndcnt equation of state (e.g., Cooray 1999; 
Garnavich ct al. 1998; Waga & Miceli 1998). Since the 
only difference between this study and previous ones is 
that we now allow w to vary with z, formalisms presented 
in previous papers should also hold except for the fact that 
w is now redshift dependent. We refer the reader to previ- 
ous work for detailed formulae and calculational methods. 
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Fig. 1. Current constraints on w(z) « w + z{dw / dz) . 
The dot-dashed lines are the current type la supernovae 
(Ricss et al. 1998; Perlmutter et al. 1998) constraints at 
the 95% confidence level, while the upper limits from grav- 
itational lensing statistics (Cooray 1999) is shown with 
solid lines. The age of the Universe as a function of Hata 
is show by dotted lines. A conservative lower limit on this 
parameter is 0.8 (dashed line) when age of the Universe 
is - 14 to 15 Gyr and H ~ 65 km s" 1 Mpc" 1 . The 
shaded region defines the parameter region allowed by cur- 
rent data. 

In Fig. 1 we summarize current constraints on w(z) as 
given in Equation 1 . Here we have assumed a flat Universe 
with Qm = 1/3. Evidence for such low matter density, in- 
dependent of the nature of the additional energy compo- 
nent, comes primarily from the mass density within galaxy 
clusters (e.g., Evrard 1997). As shown in Fig 1, there is 
a wide range of possibilities for w(z), and the allowed pa- 
rameter space is consistent with the degeneracies discussed 
in literature (e.g., Zlatev et al. 1998). Even though the 
w ~ —1/3 model has been ruled out by combined type 
la supernovae, gravitational lensing and globular cluster 

ihttp:/ /www. sdss.org 



ages, we now note that a model in which w„ ~ —1/3 but 
(dw/dz) ~ —0.9 is fully consistent with the current ob- 
servational data. 

In order to test whether one can constrain w(z) better 
than the current data, we increased the Type la supernova 
sample between redshifts of 0.1 to 1; however, the degen- 
eracy between w and (dw/dz) did not change apprecia- 
bly. Increasing the upper redshift of supernova samples 
decreased the degeneracies; thus, cosmological probes at 
high redshifts are needed to properly distinguish rcdshift- 
dependent w(z) from a constant w model. A probe to a 
much higher redshift is provided by the CMB anisotropy 
data; however, as pointed out in Hutcrer & Turner (1998) 
and Huey et al. (1998), CMB anisotropy is not a strong 
probe of w(z). This is due to the fact that w(z) affects 
mostly the lower multipoles, which cannot be measured 
precisely due to cosmic variance. Also, supernovae and 
galaxy clusters are unlikely to constrain w(z) in the near 
future given that current observational programs are not 
likely to recover them at high redshifts (z J> 2). 
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Fig. 2. Expected number of multiply imaged quasars 
with image separations between 1 and 6 arcsecs in the 
SDSS data. We have only counted lensed quasars with at 
least two images greater than a magnitude limit of 21. The 
expected number is a constant along the — lines. 
For Q M — 0.3 and tt\ — 0.7, we expect ~ 2000 lensed 
sources to be detected within the SDSS. 

Thus, an alternative probe to high redshifts is needed. 
In the next section we consider the possibility of using 
gravitational lensing statistics, in particular the redshift 
distribution of strongly lensed sources, as a probe of w(z). 
Such statistics in principle probe the volume of the Uni- 
verse out to a redshift of ~ 5. In the past, lensing statis- 
tics were hampered by the lack of large samples of lensed 
sources with a well known selection function and their 
redshift distributions, which are all needed to constrain 
w(z). An exciting possibility is now provided by the up- 
coming high-quality data from the Sloan Digital Sky Sur- 
vey 1 (SDSS; Gunn & Knapp 1993) which is going to image 
7r steradians of the sky down to a 1-a magnitude limit of 
~ 23. Since no high redshift (z ^ 2) cosmological probes 
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yet exist, gravitational lcnsing statistics may be the prime 
candidate to study w(z). 

3. GRAVITATIONAL LENSING STATISTICS 

In order to calculate the expected number of lensed 
quasars, in particular considering the SDSS, we extend 
previous calculations in Wallington & Narayan (1993; also, 
Dodelson & MacMinn 1997) and Cheng & Krauss (1998). 
Our calculation follows that of Cooray et al. (1999) in 
which we calculated the number of lensed galaxies in the 
Hubble Deep Field. We follow the magnification bias (e.g., 
Kochanek 1991) calculation in Cheng & Krauss (1998). 
We calculate the expected number of lensed sources as a 
function of the magnitude limit of the SDSS and select 
sources with image separations between 1 and 6 arcsecs. 
This range is selected based on image resolution and source 
confusion limits. Our prediction, shown in Fig. 2 as a func- 
tion of f^M and J^a, is based on current determination of 
the quasar luminosity function, which is likely to be up- 
dated once adequate quasars statistics are available from 
the SDSS. 

As shown, there are about 2000 lensed quasars down to 
a limiting magnitude of 21 that could in principle be de- 
tected from the SDSS data if Am = 0.3 and f^A = 0.7. 
This number drops to about ~ 600 when Qm = 1-0; how- 
ever, such a cosmological model is already ruled out by 
current data. In this calculation, we have ignored effects 
due to extinction and dust; this is an issue with no consen- 
sus among various studies (Malhotra et al. 1997; Falco et 
al. 1998). It is likely that a combined analysis of statistics 
from ongoing lensed radio sources and optical searches may 
increase knowledge on such systematic effects. The SDSS 
lensed quasars lie in rcdshifts out to z ~ 5. Detection of 
such high rcdshift lensed sources is likely to be aided by 
the 5 color imaging data, including the z-filter. This pos- 
sibility has already been demonstrated by the detection of 
some of the highest redshift quasars known today using 
the SDSS first year test images (Fan et al. 1998). 

4. LENSING CONSTRAINTS ON W(Z): PROSPECTS 

In order to estimate the accuracy to which one can con- 
strain w(z) based on gravitational lensing statistics, we 
take a Fisher matrix approach with parameters w and 
(dw/dz) , and assume a flat universe. As stated in liter- 
ature (e.g., Tegmark et al. 1997), Fisher matrix analysis 
allows one to estimate the best statistical errors on param- 
eters calculated from a given data set. The Fisher matrix 
F is given by: 



TP. . 



d 2 InL 

dpidpj 
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where L is the likelihood of observing data set x given the 
parameters p\ . . . p n . 

We bin the observations (number of lensed quasars) in 
rcdshift bins of width Az out to a maximum redshift z maK . 
Since the selection function for quasar discovery process is 
still unknown, we adopt a Poisson likelihood function at 
each redshift bin Az centered at redshift z. The expected 
number of lensed sources, N exp , in each redshift bin takes 
into account the magnitude limit, range of allowed image 
separations as well as the magnification bias. With these 
approximations, we can now write the Fisher matrix as: 



^ = E 



i 



dN exp (z,Az) dN exp (z,Az) 



N exp (z,Az) 



dpi 



(3) 



This form for the Fisher matrix is identical to the one de- 
rived from a Gaussian distribution when the uncertainty 
(a) of the distribution is taken to be equal to the shot-noise 
term ^/ N cxp (z, Az). 
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Fig. 3. Constraints (2-cr) on w and (dw/dz) using 
lensed source redshift distribution expected from the SDSS 
for a fiducial cosmological model of w a = —1, (dw/dz) = 
0, and marginalised over Qm (^m = 0.3 ±0.1 and = 
1-Om). We show the expected 2a errors when the redshift 
distribution of lensed sources is known with an accuracy of 
0.1 and 0.3 while the maximum rcdshift probed by lcnsing 
statistics is 3 and 5. 

In Fig. 3 we show the expected (2-er) uncertainties in 
w and (dw/dz) . Here we considered a flat universe with 
three parameters, ^m, w , and (dw/dz) . We marginalised 
over Aw allowing CIm = 0.3 ± 0.1 (1-er) and considered a 
fiducial model with w = —1 and (dw/dz)„ = 0. The three 
curves in Fig. 3 show variation of the constraint region with 
the redshift bin width Az (which is roughly equal to the 
uncertainty in the redshift determination) and with the 
maximum redshift of detected quasars z max . Photometric 
rcdshifts now allow redshift determinations with an accu- 
racy of order 0.1 68% of the time and of order 0.3 100% of 
the time (e.g., Hogg et al. 1998). It is apparent that the 
size of the constraint region decreases significantly with 
increasing z max and decreasing Az. Note that one can 
quite safely assume that z max « 5 for quasars in a survey 
such as the SDSS. Additionally, these calculations assume 
a limiting magnitude of 21, much lower than the expected 
1-er limiting magnitude of 23 for the SDSS data. 
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Fig. 4. Constraints (2-er) on w and (dw/dz)„ using 
lensed source redshift distribution expected from the SDSS 
for four fiducial cosmological models when redshift of 
lensed sources is known with an accuracy of 0.3 out to 
a redshift of 5. In all cases we assumed a flat universe and 
marginalised over Q,m (^m = 0.3 ± 0.1). 

Fig. 4 shows constraints in the w - (dw/dz) plane 
for four fiducial models. Models shown are the cosmo- 
logical constant (w(z) = —1), non-Abelian cosmic strings 
(w(z) = —1/3; Spergel & Pen 1997) and two quintessence 
models exhibiting a variation in w at small z (w(z) = 
—0.5 + 0.1 z and w(z) = —0.5 — 0.05 z). In all cases 
we assumed a flat universe and marginalised over Qm 
(Qm = 0.3 ± 0.1). We also assumed that redshifts of 
lensed objects are determined with an accuracy Az = 0.3 
and that we have data out to redshift of z max = 5. This 
figure shows that the strength of the constraints depends 
strongly on the fiducial model. In fact, we found that 
fiducial models for which w + z(dw / dz) ~ — 1 (for z of 



order unity) give weaker constraints. This result is not 
surprising and can be understood by simply using the fact 
that number of lensed objects out to a redshift of z is 
roughly proportional to the volume of the universe V(z). 
Since dV/dpi (where pi is either w or (dw/dz) ) is pro- 
portional to (1 + z Y+w a +z(dw/dz) ^ we gee that the ex _ 

pected number of lensed quasars varies slowly with pi if 
w a + z(dw/dz) ~ —1. In that case, our constraint region 
will be relatively large. 

5. SUMMARY & CONCLUSIONS 

In this Letter, we considered the possibility of constrain- 
ing quintessence models that have been suggested to ex- 
plain the missing energy density of the Universe. We sug- 
gested gravitational lensing statistics, which can be used as 
a probe of the equation of state of the missing component, 
w(z). An exciting possibility to obtain an adequate sample 
of lensed quasars and their redshifts comes from the Sloan 
Digital Sky Survey. Writing w(z) w w a + z(dw/dz) , we 
studied the expected accuracy to which equation of state 
today w and its rate of change (dw/dz) can simulta- 
neously be constrained. Adopting some conservative as- 
sumptions about the quality of the data from SDSS and 
assuming a flat universe with VLm = 0.3 ± 0.1, we con- 
clude that tight constraints on these two parameters can 
indeed be obtained. The strength of the constraints de- 
pends not only on the quality of the lensing data from the 
SDSS, but also on the fiducial model (true values of w 
and (dw / dz) a ). In particular, fiducial models for which 
w a + z{dw / dz) ~ —1 (for z of order unity) give weaker 
constraints on w and {dw/dz) a . 
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